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We report giant Goos-Hanchen shifts [F. Goos and H. Hanchen, Ann. Phys. 436, 333 (1947)] 
for electron beams tunnehng through graphene double barrier structures. We find that inside the 
transmission gap for the single barrier, the shift displays sharp peaks with magnitudes up to the 
order of electron beam width and rather small fuU-widths-at-half-maximum, which may be utilized 
to design valley and spin beam splitters with wide tunability and high energy resolution. We 
attribute the giant shifts to quasibound states in the structures. Moreover, an induced energy gap 
in the dispersion can increase the tunability and resolution of the splitters. 
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In optics, it is well-known that a light beam under- 
goes a lateral shift when it is totally reflected from a 
dielectric interface ii This phenomena is referred to as 
the Goos-Hanchen (GH) shif1j2i and can be theoretically 
explained based on the reshaping of the wave packet.— 
Analogies of the GH shift have been widely considered 
in various fields, including acoustics;^ electronics?^^— rel- 
ativistic corrections^ atomic optics?^ and neutronics 
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Recently, the analogy for the massless electron has 
become of interes1>ii"— since discovery of graphene, a 
monolayer of sp^ bonded carbon atomSf^i It has been 
reported that the GH shift plays an important role in 
the group velocity of quasiparticles along interfaces of 
graphene p-n junctions ; ^^'^^ whereby a twofold degen- 
eracy on top of the usual spin and valley degenera- 
cies can be introducedii and coherent buffers and mem- 
ories can be realized in graphene p-n-p waveguidesi^ 
The valley- dependent GH shif t^^'^^ based on strained 
graphene has been very recently utilized by Zhai et al. 
to design a valley beam splitterJ^ We note that, to ef- 
fectively realize the proposed splitter, the difference of 
GH shifts for valley K and K' should be larger than 
the longitudinal width of the electron beam defined as 
Wy = Wb cos~^ a, where Wb and a are the waist width 
and incident angle of the electron beam, respectively. 
This condition also ensures the validity of the stationary- 
phase approximatioufi^ which is widely used in GH shift 
calculationiir— especially for resonant transmission.^- 
Considering a typical beam divergence {Sa = Xp/T^Wb 
with \f the Fermi wave-length) of l°-0.1°r^ Wb is 180A_f- 
ISOOAf- In a n-p-n single barrier structure (SBS)r^ 
this condition may not be met without properly se- 
lecting the structural parameters. On the other hand, 
the spin-dependent GH shift was also proposed to spa- 
tially split spin beams, based on a SBS formed by a 
local magnetic field^S and an electrostatic potential in 
a two-dimensional electron gas (2DEG)ii However, the 
displacements of the spin-dependent GH shifts are also 
found to be insufficient<i Therefore, finding a suitable 
structure to effectively implement these ideas is one of 
the main issues in such an exciting field. 

In this letter, we report a giant GH shift of elec- 
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FIG. 1. (color onhne) (a) Sectional and (b) top schematic di- 
agrams for a particle quantum tunneling through a graphene 
DBS, with the two barriers of width toi(2), height Ui{2), and 
distance d between them. I = w\ + W2 + d is the total width 
of the structure. In (a), the dashed lines show smooth electric 
potentials with distributions of error functions and transition 
regions' widths of 0.1«;i(2). In (b), the upper (red solid) and 
lower (blue dashed) components have their locus at ±5/2, 
(7r =F '5/2, and at ± 5/2 for the incident, reflected and trans- 
mitted beams, respectively. A detector placed in a proper 
position of the outgoing region can detect the giant GH shift. 



tron beams tunneling through a double barrier structure 
(DBS) in graphene. Wc find that inside the transmission 
gap (TG) for the constituted single barrier— the GH shift 
displays sharp peaks which are absent in the SBS cases 
and were attributed to the quasibound states formed in 
the DBS. Remarkably, we find that the magnitudes of 
the peaks are much larger than the maximum magnitude 
in the corresponding SBS and can easily achieve the or- 
der of Wy. Together with the rather small full-widths at 
half-maximum (FWHM) of the peaks, this giant GH shift 
can be used for designing valley and spin beam splitters 
with wide tunability and high energy resolution. We also 
investigate the effects of the structural asymmetry of the 
DBS and the induced energy gap in the linear dispersion. 
The results show that the former suppresses the GH shift 
while the latter enhances it, which increases the control- 
lability of devices based on the GH shift. 

To investigate the GH shift in graphene DBSs, we 
consider an incident beam well collimated around some 
transverse wave vector q G {—qm,q_m) {<lm = E/hvp) 
[corresponding to a = arcsin(g/(7„i)], quantum tunneling 
through a graphene DBS as sketched in Fig. 1. The 
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Dirac equation can be written as: 

(ufcr •p + Act.. +C/)* = (1) 

where vp sa lO^m/s is the Fermi velocity, cr = {ax, <Jy) is 
the pseudospin operator given by Pauh's matrices, p = 
{Px,Py)^ is the momentum operator, A = mvp{Aso) 
is the energy gap owing to the sublattice symmetry 
breakingSS (the spin-orbit interaction^^) , and U = C/i(2) 
in the barriers while U = elsewhere. The wave packets 
of the incident and associated reflected and transmitted 
beams at the two terminals of the DBS can be given by 

(2a) 

(2b) 

= ^ ...(.)/(.-,)e-(^)--(;^,„,,/,). 

(2c) 

Here each plane wave of the spinor form is a solution 
of Eq. (1), and a basis is used to ensure that the 
product of the upper and lower components is realj^ 
f{q — q) is the angular spectral distribution which can 
be assumed to be of Gaussian profile, Wye~^y'^'^~^^ 1"^ . 
k{q) = \Jq^ — (A//i?;i?)2 — is the longitudinal wave 
vector and a'((?) = ^m~^{hvpql \J E'^ — A^) is the phase 
angle. Due to the induced energy gap A, a' is usu- 
ally different from the incident angle and a factor A = 
\/ — A'^ / [E + A) appears in the lower component of 
the spinors. For the reflected and transmitted beams, 
r{q) = |r(g)|e''>'-(«) and i(g) = |t(g)|e''^'(9) are the reflec- 
tion and transmission coefficients respectively, which can 
be determined by the continuities of each component and 
calculated by the standard transfer-matrix method)^ 

From the stationary-phase approximatioufiS the van- 
ishing of the gradient of the total phase in direction gives 
the locus of the steady wave packet peak. Due to the 
spinor nature of graphene, the loci of the two components 
in a beam are found to be different i^i For the incident 
and transmitted beam, the upper component exceeds the 
lower one with a distance of (5 ((5 = 1/A:) , while for the re- 
flected beam the inverse becomes the case [see. Fig. 1(b)]. 
The deviation between the loci of the reflected (trans- 
mitted) and the incident beams thus gives corresponding 
GH shift, ~ ~d(j)r/dq\g =F S and = —d4>t/dq\q, 
where ipt = 4>t + kl. Note the GH shift in reflection is 
component dependent while the shift in transmission is 
not [see. Fig. 1(b)]. We'd like to use ct^ = —{d(l)r/dq)\q to 
describe the average shift in reflection. 

Fig. 2(a) shows the calculated GH shift for a symmet- 
ric graphene DBS as a function of at a fixed a = 10°. 
In this calculation, A = 0. To make sure the electron 
density of states coincides with a true system and the 
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FIG. 2. (color online) (a) The GH shift (solid) in a symmetric 
graphene DBS [C/ = f/i = U2 =62 meV, wi = W2 =100 nm, 
d =50 nm] as a function of the incident energy at q = 10°. 
The insert shows the GH shift in a SBS with the same barrier 
height and width. The dashed line indicate the semi-classical 
shift predicted by the Snell's law, which has no definition in 
the TG. (b) The dependence of the shift sharp peak(s) on 
the incident angle, a = 10°, 15°, and 20°. The positions and 
heights of the four giant GH shift peaks are marked. 



system stays in a ballistic regime, typical values are used: 
U = 62 meV; w — 100 nm. and d =50 nm. These 
parameters also ensure the legitimacy of the stationary- 
phase approximationii^iiSi^^ One can see from Fig. 2(a) 
that outside the TG, the GH shift shows the same trend 
as the shift predicted by geometric optics using Snell's 
lawM^ Moreover, the GH shift oscillates around the 
Snell shift as it is enhanced (suppressed) at kiWi = Ntt 
[k,Wi = {N - l/2)7r] with N = 1,2,3,... where reso- 
nant (antiresonant) tunneling happens. Note, the posi- 
tive peak nearest the TG (denoting as Psbs) which has 
the maximum magnitude for tunneling through a SBS 
cannot be enhanced to the order of the longitude beam 
width without properly selecting the structural parame- 
ters. 

Remarkably, a significantly sharp peak with a magni- 
tude comparable to Psbs appears inside the TG [around 
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FIG. 3. (color online) The phase shift of the transmitted beam 
vs. the incident angle for beams with different E/U shown in 
the figure. (a)/(b) for a rectangular/smooth DBS. 



E/U = 1 in Fig. 2(a)]. This is absent in the case of elec- 
tron beams tunneling through graphenc SBSs [see the 
inset of Fig. 2(a)]. One can clearly see from Fig. 2(b) 
that with increasing incident angle the peak value of the 
GH shift dramatically increases. At a = 20°, the peak 
value reaches ^ 1700Af (~ 16000 nm) which is about 
ten times of the corresponding Psbs and is in the order 
of the longitudinal beam width. Therefore, the obtained 
results suggest that the valley splitter based on this struc- 
ture can be realized with much looser conditions, since 
giant GH shift for valley K{K') can be obtained inside 
the TG while the GH shift for valley K'{K) retains the 
order of Xp- In addition, due to the rather small FWHM 
the splitter will possess a much higher energy or wave 
vector resolution. 

It is surprising that the GH shift in a graphene DBS 
possesses such giant magnitudes, since the shift along 
a single interface or through a SBS (inside the TG) is 
of the order of Fermi wave-length Afi^"— Here, we at- 
tribute it to the quasibound states in the DBS, which are 
formed by the evanescent waves in the two barriers. It 
is well-known that plane waves with different q generally 
present different phase shifts, leading to the reshaping of 
the wave packet and thus the GH shift of the electron 



FIG. 4. (color online) (a) GH shifts in transmission and re- 
flection for an asymmetric DBS with U =62 meV, Ui = 0.8f/, 
U2 = 1-2U, wi ^ W2 ^ w =100 nm, and d =50 nm. (b) GH 
shifts in transmission for different cases of structural asym- 
metries with the parameters (wi / w , W2 / w) and {Ui/U, U2/U) 
in the sequence given in the figure. For all cases, a = 10° . 



beam. In our case, when the center plane wave is aligned 
to the quasibound state, multiple interferences will arise 
through the quasistanding waves between the two barri- 
ers. This will lead to a remarkable difference in the phase 
shifts between the center plane wave and adjacent ones 
[sec Fig. 3(a)]. Accordingly, a giant GH shift for such 
an electron beam is present inside the TG. One can also 
see clearly from Fig. 3(a) that with the increasing E/U 
of the beam, the incident angle which shows the biggest 
slope of ~d4>r/da also increases, which is consistent with 
Fig. 2(b). Thus, the observed giant GH shift inside the 
TG is due to the quasibound states formed in graphene 
DBSs. 

One may wonder how the smoothness of realistic bar- 
riers will affect the magnitude of the GH shift. To 
answer this problem, we consider a realistic potential 
which varies smoothly on the scale of the graphene lat- 
tice constant and adopt a typical potential profile^ of 
U{x) = 0.5C/[erf(2x/Lb-2)-herf(2(w-a;)/i6-2)], where 
erf(x) is the error function and the width of the transi- 
tion region is set as Lb = O.lw (see Fig. 1(a)). We cal- 
culated the transmission coefficient, phase shift, and also 
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FIG. 5. (color online) The GH shift in transmission through 
a graphene DBS with wi = W2 =50 nm, d =100 nm and 
U = Ui=U2 =50 meV. The induced gap is A = 0, 5, 10, 15 
meV as indicated in the figure. Insert: the GH shifts for the 
refiected and transmitted beams for A = 5 meV in the same 
structure. For all cases, ct = 20°. 

the GH shift. Similar results have been achieved. Fig. 
3(b) shows the phase shifts as a function of incident angle 
for a smooth DBS. Compared with the case of a rectan- 
gular DBS with the same a = 15°, the giant GH shift for 
a smooth DBS appears at a lower incident energy (from 
E/U^lMl to 0.946) and the magnitude of the giant GH 
peaks decreases from 245. 4Ai? to 67.4Af. To achieve the 
GH shift in the order of the longitude beam width, one 
can increase the incident angle (about 23°). Thus the 
smoothness of the potential barriers will not restrict the 
use of the giant GH shift. For convenience, we will adopt 
the rectangular model in the following discussions. 

We now study the effect of the structural asymmetry of 
the DBS which may be introduced during sample prepa- 
ration and gate control, i.e., the two barriers have differ- 
ent widths and/or heights. Fig. 4 shows the GH shifts in 
transmission and reflection for asymmetric DBSs. Inter- 
estingly, one can see from Fig. 4(a) that the GH shift in 
reflection crucially depends on which terminal the beam 
is being reflected from, where tXr and at (cr^ ^-nd "'t) mean 
the GH shift in reflection and transmission for electron 
beams incident on the Ui (U2) terminal (as shown in Fig. 
1). We would like to point out that Ur and cr^ arc found 
to follow the relation of ar + cr'^ = 2(Tt which holds for 
any asymmetric DBSs and can be understood by the rela- 
tion between the reflection and transmission coefficients. 
Through the scattering matrix;^ we can get t' = t and 
r' ~ ^tr* /t*, which imply = (pt and (pr + (j)'^ — 20t, 
respectively, r (r') and t {t') are the reflection and trans- 
mission coefficients of the beam reflected from or tunnel- 
ing through the Ui (U2) terminal, where the term kl has 
been contained in t and t' . Using the definition of GH 
shift, we get at = a't and ar + a'^ = 2a t- Note for a sym- 
metrical DBS, 04 = (f>r {4>'t ^ 4''r) thus 0^ = (j)r, which 



indicates that the two GH shifts in reflection will become 
identical in this case. 

We now take at {a[) as a target to evaluate the ef- 
fect due to the structural asymmetries. As shown in 
Fig. 4(b), the effect of structural asymmetry always sup- 
presses the GH shift inside the TG. For the case of barri- 
ers with the same height but different widths, the struc- 
tural asymmetry only decreases the magnitude of the GH 
shift inside the TG while not moving the peak position. 
For the case of barriers of different heights, the effect of 
structural asymmetry not only decreases the magnitude 
of the GH shift inside the TG but also makes the peak 
position shift. Moreover, the FWHM of the peak of the 
GH shift increases. 

Fig. 5 shows the GH shift for the transmitted beam 
in a graphene DBS with different induced energy gaps in 
the linear dispersion. As is seen, when there is a nonzero 
gap, the GH shift has no definition for cos^ a < 
as the electron cannot propagate freely even in the non- 
modulated regions. Meanwhile, the GH shifts for the 
reflected and transmitted beams differ (see insert in Fig. 
5). With increasing induced energy gap, the peak posi- 
tions of the GH shift inside the TG move to higher en- 
ergies. Moreover, the magnitudes of the peaks increase 
and the FWHMs become even smaller, which increases 
the tunability and energy resolution of the valley or spin 
splitter device. The underlying physics is that, the pres- 
ence of energy gap increases the modulus of the longitu- 
dinal wave vector (/t), which makes a stronger multiple 
interference effect and thus a bigger GH shift at a little 
higher energy. 

In summary, we have theoretically calculated the GH 
shifts of reflected and transmitted electron beams in a 
graphene DBS. Interestingly, we found that the GH shift 
displays sharp peaks inside the TG for the constituted 
SB, which are absent in the graphene SBS cases and 
can be attributed to the quasibound states formed in the 
DBS. The reported giant GH shift can be detected in the 
transmitted beam by placing a detector in the outgoing 
region and far away from the incident position (see Fig. 
1(b)). In the TG, the detector will collect no electrons 
unless the incident energy is aligned to the quasibound 
states. The results obtained in this work suggest the fea- 
sibility of making valley splitter based on graphene DBSs. 
The tunability and energy resolution of the splitter can 
be further increased by an induced energy gap in the lin- 
ear dispersion. By the spin-dependent giant GH shift, 
a spin beam splitter restricted by the small magnitude 
of the GH shiftii may be realized in 2DEG or graphene 
based DBSs now. 
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